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Abstract. In the framework of the dwarf nova thermal- 
viscous disc instability model, we investigate the combined 
effects on the predicted dwarf nova lightcurves of irradi- 
ating the accretion disc and the secondary star and of 
evaporating the inner parts of the disc. We assume the 
standard values of viscosity. We confirm the suggestion 
by Warner ( |1998| ) that the large variety of observed out- 
bursts' behaviour may result from the interplay of these 
three effects. We are able to reproduce light curves rem- 
iniscent of those of systems such as RZ LMi or EG Cnc. 
We can obtain long lasting outbursts, very similar to su- 
peroutbursts, without assuming the presence of a tidal 
instability. 

Key words: accretion, accretion discs - instabilities - 
(Stars:) novae, cataclysmic variables - (stars:) binaries : 
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1. Introduction 

Dwarf novae (DN) are cataclysmic variable binary systems 
which, every few weeks, exhibit 4-6 mag outbursts, which 
last for a few days (see e.g. Warner 1995a). In several 
subclasses of DN both the outburst durations and recur- 
rence times can be very different from the values quoted 
above. It is generally accepted that DN outbursts are due 
to a "thermal-viscous" instability. This instability occurs 
in the accretion disc in which the viscosity is given by 



the a-prescription (Shakura & Sunyaev 1973), at temper- 
atures close to 8000 K. Hydrogen is then partially ionized 
and o paciti es are a steep function of temperat ure (s ee Can- 



1993 for a review and Hameury et al. 199g for the 



most recent version of the model). Modeling dwarf nova 
lightcurves requires a varying Shakura-Sunyaev parameter 
a (Smak |1984|): it must be of the order of 0.1 - 0.2 (0.2 



according to Smak 1999b) in outburst and of the order of 
0.01 in quiescence, when the temperature is below the hy- 
drogen ionization temperature. One could therefore expect 
that the disc instability model (DIM) might offer useful 
constraints on mechanisms which generate accretion disc 
viscosity. This assumes, of course, a successful application 
of the model to the observed DN outburst cycles. 

However, despite its success in explaining the overall 
characteristics of DNs, the DIM in its standard version 
faces several serious difficulties when one tries to account 
for the detailed properties of dwarf nova outbursts. Some 
of these difhculties are the result of an incomplete version 
of the DIM. For example it was believed that a trunca- 
tion of the inner parts of the disc is necessary to explain 
the long delay between the rise of optical light and that 
of UV an d EUV in systems such as SS Cyg. As shown 
by Smak ( 1998 ), however, when correct outer boundary 
condition are assumed, the standard DIM reproduces the 
observed delays. On the other hand, observed quiescent 
X-ray fluxes far exceed the predictions of the model and 
seem to require an inner 'hole' in the disk. Such a hole 
can either be due to evaporation of the disc close to the 



white dwarf (Meyer & Meyer-Hofmeister 1994), or to the 
presence of a magnetic field strong enough to disrupt the 
disc. In addition, systems such as WZ Sge, which have 
long recurrence time and large amplitude, long outbursts 
require very low values of a (a < 10^*) if interpreted in 
the framework of the standard DIM (Smak [1993 , Osaki 
1995a, Meyer-Hofmeister et al. 1998| ). These values, much 
lower than those of other DNs at similar orbital periods, 
are, however, left unexplained. On the other hand, WZ 
Sge systems can be explained with standard values of a, 
provided that the disc is truncated as in other systems so 
that it is either stable or marginally unstable (Lasota et al. 



1995, Warner et al. 1996) and the mass transfer from the 
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secondary is significantly increased during the outburst 
under the influence of illuminatio n by radiation from the 
accreting matter (Hameury et al. 1997 ). 
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SU UMa systems are a subclass of dwarf novae which 
occasionally show long outbursts during which a light- 
curve modulation (superhump) is observed at a period 
slightly longer than the orbital period; these superout- 
bursts are, usually, separated by several normal outbursts. 
The superhump is due to a 3:1 resonance in the disc which 
causes the disc to become eccentric and to precess (White- 
hurst 1988). Osaki ( |1989D proposed that the related tidal 
instability is also responsible for the long duration and 
large amplitude of superoutbursts (see Osaki 1996 for a re- 
view of the thermal-tidal instability model) . In his model, 
the tidal torques which remove angular momentum from 
the outer parts of the disc are increased by approximately 
an order of magnitude when the disc reaches the 3:1 reso- 
nance radius (typically 0.46a, where a is the orbital sepa- 
ration) until the disc outer radius has shrunk to typically 
0.35a, i.e. by about 30%. This model accounts for many 
of the properties of SU UMa systems; it has, however, dif- 
ficulties in explaining systems with very short superout- 
burst cycles such as RZ LMi, for which one must assume 
that the tidal instability stops when the disc has shrunk 



by less than 10% (Osaki 1995b), i.e. by much less than 
assumed in the standard case. 

Finally, it must be noted that a number of systems ex- 
hibit very bizarre lightcurves: we have already mentioned 
the case of RZ LMi, which has similarities, with systems 
such as ER UMa, VI 159 Ori and DI UMa. The case of the 
December 1996 superoutburst of EG Cnc which was fol- 
lowed by 6 closely spaced normal outbursts in 1996 has not 



been reproduced by simulations, except by Osaki (1998) 
who assumed that the viscosity parameter ac in the cold 
state was increased to 0.1, almost the value in the hot 
state, for 70 days after the superoutburst, and then re- 
turned to its quiescent value ac = 0.001. It is finally worth 
mentioning that the prototypical classical dwarf nova U 
Gem exhibited an unusually long outburst in 1985, last- 
ing 45 days, with a shape similar to superoutbursts, but 
without superhumps (Mattel et al. 1987 ). Since in this 
system the radius of the 3:1 resonance is larger than the 
size of the primary's Roche lobe, one can conclude that 
while superhumps can be attributed to the 3:1 resonance, 
the tidal instability is obviously not the sole cause of very 
long outbursts. 

The inability of numerical models to reproduce the 
large variety of observed light curves may indicate that 
additional physical effects should be added to the DIM. 
One such effect is the tidal instability. Another impor- 
tant class of effects is the illumination of the disc and the 
secondary star. These effects are usually not included in 



simulations, and it was suggested by Warner (1995c) that 
irradiation of the secondary star, that gives rise to high 
and low states of mass transfer, and of the inner disc, that 
drastically affects the disc instability, account for the light 



secondary to irradiation. This suggestion is supported by 
observations since there is evidence that the mass trans- 
fer rate from the secondary star increases during outbursts 
(Smak |l995| ) in dwarf novae such as Z Cha or U Gem, most 
probably as a result of illumination of the secondary. 



Recently Smak ( |l999a| ) argued that properties of out- 
burst cycles of "standard" U Gem-type dwarf novae, 
which cannot be reproduced by the DIM (e.g. the same 
maximum brightness of narrow and wide outbursts), are 
well explained if all outbursts are associated with some 
mass-transfer enhancement due to the secondary's irradi- 
ation. (Long outbursts would be due to important mass 
transfer enhancements, making the mass transfer rate 
larger than the critical value for stability). If this is the 
case the 'pure' DIM would find no application in the real 
world. 

It also appears that illumination of the disc itself by the 
hot white dwarf has strong effects on the stability prop- 
erties of the disc as soon as the white dwarf temperature 
exceeds 15,000 K (King |1997| , Hameury at al. |1999| ). 

In this paper, we investigate the influence of the com- 
bined effects of illumination of the disc, of the secondary 
star and of evaporation of the inner parts of the disc, on 
the predictions of the DIM in which standard values of 
viscosity are assumed. Our free parameters are the white 
dwarf mass, the quiescent white dwarf temperature, the 
mass transfer rate in the absence of illumination, and two 
parameters describing in a crude manner evaporation ef- 
fects and the influence of illumination of mass transfer 
from the secondary. We show that a large variety of light 
curves are predicted by the models, many of which have 
an observational counterpart. In section 2, we show that 
the long 1985 outburst of U Gem requires enhanced mass 
transfer during the outburst. In section 3, we describe the 
model and our assumptions; in section 4, we give our re- 
sults and compare them with lightcurves of observed sys- 
tems, and we discuss briefly possible extensions of this 
work. 



2. The long 1985 outburst of U Gem 

U Gem is a prototypical dwarf nova that undergoes out- 
bursts that last for 7-14 days, with an average recurrence 
time of 120 days (Szkody & Mattel |1984| ). The orbital pe- 
riod is 4.25 hr, and the primary and secondary masses are 
respectively 1.1 and 0.5 M© (Ritter & Kolb|1998D; for such 



parameters, the outer disc radius, defined as the radius of 
the last non-intersecting orbits, is 4.15 10^° cm. 

The 1985 outburst that lasted for 45 days is there- 
fore exceptional; long outbursts are observed in other sys- 
tems (SS Cyg for example), and are a natural outcome 
of models in which the outer disc radius is kept constant 



curves of VY Scl stars. He also suggested (Warner 1998) 



(see e.g. Hameury at al. 1998), which may happen when 



that the wide spectrum of superoutburst behaviours is 
generated by the interplay of reactions of the disc and the 



the outer disc radius reaches the tidal truncation radius. 
However, one does not normally obtain extremely long 
outbursts; the 1985 outburst of U Gem was in fact so long 
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that, more mass was accreted during this outburst than 
what was contained in the pre-outburst disc. The maxi- 
mum mass Md,max of a quiescent disc during quiescence is 
the integral of the maximum surface density Smax on the 
cool branch; using the fits given by Hameury et al. ( 1998D , 
one gets 



2 X 10^"a 



20„ -0.83 j,f- 



?'out 

IQiOcm 



3.14 



(1) 



where etc is the Shakura-Sunyaev viscosity parameter on 
the cool branch, Afi is the primary mass in solar units, 
and Tout is the outer disc radius. 

On the other hand, during a long outburst, the whole 
disc is entirely on the hot branch, and the local mass trans- 
fer rate in the outer regions of a disc must be large enough 
to prevent a cooling wave formation; using again the ana- 
lytical fits of Hameury al al. (1998), this gives 



> 



^ VlQiOcm/ 



2.67 



(2) 



Here, Afout is the local mass transfer rate in the outer 
parts of the disc, which is close to the mass accretion rate 
onto the white dwarf if the whole disc sits for a long time 
on the hot branch. 

The maximum duration of such an outburst is thus 
^max = Afd,max/A/. For the parameters appropriate for U 
Gem, one gets 



in 



26 



VO.Ol/ ^ V 



4.1 10i"cm 



0.47 



(3) 



As Qfc is larger than 0.01 for this prototypical dwarf nova 
(Livio & Spruit 1991 for example find that ac must be 
equal to 0.044 to account for the timing properties of U 
Gem), trminax Can never be as high as 45 days. This means 
that the total amount of mass accreted during this long 
outburst is larger than the mass of the disc in quiescence; 
this is possible only if the mass transfer rate from the sec- 
ondary has increased to a value close to the mass accretion 
rate onto the white dwarf, i.e. is close to the critical rate for 
stable accretion. Such an increase of the mass transfer rate 
is very likely caused by the illumination of the secondary. 
This reinforces the conclusion of Smak (1999a) that long 
outbursts result from large mass-transfer enhancements. 

3. The model 

3.1. Disc irradiation 

We use here the numerical code described in Hameury et 
al. (1998). This code solves the usual mass, angular mo- 
mentum and energy conservation equations on an adap- 
tive grid, with a fully implicit scheme. This allows to re- 
solve narrow structures in the accretion disc (Menou et 
al. 1999), and avoids the Courant condition which would 



severely limit the time step. To describe disc irradiation 
we us e a version of the code d escrib ed in Dubus et al. 
( |l999|) (see also Hameury et al. 1999 ). A grid of vertical 



structures is used to determine the cooling rate of the disc 
as a function of the vertical gravity, the integrated disc 
surface density S, the central temperature Tc and the il- 
lumination temperature Tm, defined as Tm = {Fm/aY^^ 
where -Flu is the illuminating flux. In what follows, we 
use aiiot ==0.2 and acoid = 0.04, except where otherwise 
stated. We also neglect the albedo (3 of the disc; taking it 
into account introduces a multiplicative factor (1 — (3)~^^'^ 
for the white dwarf temperatures. 

It must also be stressed that the white dwarf surface 
temperature cannot be too large; this is because the in- 
trinsic (quiescent) white dwarf luminosity must be signif- 
icantly less than the accretion luminosity in outburst; for 
Ml — 0.6 Mq, the quiescent white dwarf temperature has 
to be smaller than 33,000 K if the outburst amplitude is 
to be larger than 2 magnitudes. 

3.2. Illumination of the secondary 

We are interested here in the effect of illumination of the 
secondary on its mass transfer rate on short time scales 
(days), and we do not consider any long term effects that 
may lead to cycles accounting for the observed dispersion 
of the average mass tr ansfe r rate for a given orbital period 
(McCormick & Frank [l99^ . Even on short time scales, the 
response of the sec ondar y to illumination is complex (see 
e.g. Hameury et al. 1988); we prefer to use here a simpler 
approach in which we assume a linear relation between the 
mass transfer rate from the secondary A/tr and the mass 
accretion rate onto the white dwarf Afacc, i-e. 

Mtr = max(Afo,7A/acc) (4) 

where Mq is the mass transfer rate in the absence of illumi- 
nation. This is similar to the formula used by Augusteijn et 
al. (1993) in the context of soft X-ray transients. Although 
it is an extremely crude approximation, it has, neverthe- 
less, the advantage of having only one free parameter 7. 
Its value must be in the range [0 — 1] for stability reasons. 

Such an approach obviously requires the illumination 
to have a noticeable effects on the secondary's surface lay- 
ers. As mentioned earlier, strongly irradiated companion 
stars are observed in several systems. To describe the ef- 
fects of irradiation of a Roche-lobe filling star we shall 
follow the approach of Osaki (1985) and Hameury et al. 
(|198§). 

The mass transfer rate fr om the secondary can be writ- 
ten as (Lubow & Shu |l975|) : 

A/tr = QPLlCs (5) 

where Q is the effective cross section of the mass transfer 
throat at the Lagrangian point Li, g = 1.9 X 10"T4Ph2, 
cm^ , where T4 is the surface temperature of the secondary 
and Phr the orbital period in hours; Cg is the sound speed, 
and pli the density at Li, which, in the case of an isother- 
mal atmosphere, can be expressed as: 

,{R-Rli)/H 



Pli = Poe^ 



(6) 
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where po is the density calculated at a reference level, R 
the secondary radius (defined by this reference position), 
Rli the Roche lobe radius and H the scale height. 

In quiescence, the mass transfer rate is low, of or- 
der of 10^^ — 10^^ g s~^. At short orbital periods and 
correspondingly low secondary's temperatures, both Q 
and Cs are small. If one assumed that po corresponds to 
the photospheric density, which for a low mass main se- 
quence star is of the order of 10^^ g cm^'^, one would find 
Mtr ~ ioi6g(fl-flii)/-ff g g-i. This in turn would mean 

that [R — Rli) / H is small so that Mtr varies as T^^^, and 
therefore is not very sensitive to illumination. In the ir- 
radiated case, however, the reference density po does not 
correspond to the photospheric density but to the density 
at the base of the isothermal atmosphere, which extends 
much deeper than the unilluminated photosphere when 
the atmosphere is illuminated with an irradiation temper- 
ature exceeding lO** K. For those high illumination fluxes, 
the outer layers of the star are affected on a thermal time 
scale (seconds) at least down to a point where the unper- 
turbed temperature equals the irradiation temperature. 
This point is the base of the isothermal layers in the il- 
luminated case, and from models of very low mass stars 



Tin is a simple function of the mass accretion rate onto the 
white dwarf: 



(Dorman et al. 1989| ), one gets pq ^ IQ ^ g cm with re- 



sulting very high mass transfer rates, exceeding 10 g s . 

The latter estimate is an upper limit, as it does not 
take into account the fact that a fraction of the secondary 
is shielded from irradiation by the accretion disc; par- 
tial shielding does not suppress the enhancement of mass 
transfer, since circulation at the surface of the star pre- 
vents the existence of large temperature gradients, but 
reduces it in a complex way which we are not attempting 
to describe here. 

Observations show, however, that the increase can be 
quite significant; the mass tra nsfer rate rises by a factor 2 
in U Gem and Z Cha (Smak |l995| ), whereas Vogt ( |l983|) 
finds that in VW Hyi, the bright spot luminosity increases 
by a factor ~ 15 during maximum and decline of outbursts 
close to a superoutburst, which he attributed to a corre- 
sponding increase in mass transfer under the effect of illu- 
mination. Although the evidence is not very strong, there 
are some indications that the hot spot brightening, and 
hence the mass transfer increase, is delayed with respect 
to the eruption by a day or two; this could be either the 
response time of the secondary (Smak 1995| ), or the ther- 
mal inertia of the white dwarf (only an equatorial belt is 
instantaneously heated by accretion, so irradiation of the 
secondary is delayed). 



3.3. Inner disc radius 

There is evidence that in dwarf novae accretion discs are 
truncated, as indicated by emission line profiles (Men- 
nikent & Arenas 1998), or the detec tion o f a significant 



Tin = 9.8 X l(fM^^''^M[ 



1/7 4/7 

M30 cm 



(7) 



where M15 is the mass accretion rate in units of lO^'"^ g s~^. 
Ml is the white dwarf mass and p,^Q is the magnetic mo- 
ment in units of 10^° Gcni^. The value of should be 
such as to allow ri„ = Ri in outbursts, where Ri is the 
primary radius, as most DNs do not show then coherent 
pulsations. In quiescence, however, coherent oscillations 
are observed (Patterson et al. 1998|), and for example in 
WZ Sge « 50 (Lasota et al. |1999D . 

An inner hole in the disc can be also due to evapora- 
tion. The physics of evaporation is poorly understood but 
sever al models were proposed (Meyer & Meye r- Hofmeister 



1994 , Liu et al. |l997i Kato & Nakamura |998|, Shaviv et al. 
1999 ). The evaporation rates E in the disc are, however, 



quite uncertain. Evaporation is normally accounted for by 
introducing the additional term E in the mass conserva- 
tion equation. However, because evaporation is expected 
to increase towards the accreting body, and since the local 
mass transfer rate in the disc increases sharply with radius 
during quiescence, the effects of evaporation are important 
essentially very close to the disc inner edge, and can be 
treated assuming that the disc inner radius is a function 
of the accretion rate onto the white dwarf, just as in the 
case of the formation of a magnetosphere. To first order, 
the effect of evaporation is to create a hole in quiescence 
which increases the recurrence time; what matters is thus 
the inner disc radius in quiescence and not the detailed 
way Tin varies with Mace- 

We shall therefore use equation (Q) in all cases, using 
/i30 as a free, unconstrained parameter that merely de- 
scribes the size of the hole generated in the disc by either 
the presence of as magnetic field, or by evaporation. 

4. Results 

In the following, we discuss the influence of each indi- 
vidual effect mentioned above; our reference situation is 
that of a system with a 1.0 Mq primary, whose radius is 
5 X 10^ cm; the mass transfer rate is 3 x 10^^ g s~^ , and the 
average outer disc radius is 1.8 x 10^° cm. These parame- 
ters are typical of short period dwarf novae with massive 
primaries. The light curve corresponding to the standard 
version of the DIM is given in Fig. (|l|) . For comparison we 
also show the case of a system with a 0.6 M© primary. 

Influence of the disc irradiation 

The effects of irradiation of the disc by both the hot white 
dwarf and the boundary layer has been described in de- 



tail in Hameury et al. (1999), and we summarize here the 



quiescent X-ray and UV flux (Lasota |l996 ). If this is due 
to the presence of a magnetic field, the inner disc radius 



most important results. For very hot white dwarfs (Tcff > 
20 000 K), the temperature in the innermost parts of the 
disc exceeds the hydrogen ionization temperature during 
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40 

time (days) 

Fig. 1. Predicted light curve in the standard model (no ir- 
radiation, fixed inner radius) for an average mass transfer 
rate of 3 x 10^^ g s~^, an average disc radius of 1.9 x 10^" 
cm, and a 0.6 Mq (top panel), or 1.0 Mq (bottom panel) 
primary. 
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Fig. 2. Predicted light curve when disc illumination is 
taken into account. The parameters are the same as in 
Fig. |l|. The upper panel shows the case of a 0.6 M© pri- 
mary in which the radiation from the boundary layer is 
ignored, the lower panel corresponds to a IMq primary, 
and we have taken into account the accretion luminosity. 
In both cases, the quiescent white dwarf temperature is 
30,000 K. 



quiescence; the viscosity is therefore high in these regions, 
which are thus partially depleted as first suggested by 
King (1997). The transition region between the hot in- 



ner disc and the outer, cool parts is strongly destabilized 
by irradiation, and the model predicts several small out- 
bursts between major ones. In particular, many reflares 
are expected at the end of a large outburst (see Fig. ^ . In 
certain cases, the reflares may dominated the light curve; 
this depends on whether the heating front can reach the 



outer edge of the disc or not. The reflares we obtained do 
not have the observed amplitudes but it is tempting to 
attribute the succession of several normal outbursts after 
a superoutburst in EG Cnc to this effect. Playing with pa- 
rameters would produce a result corresponding better to 
the observed lightcurve but the merit of such an exercise 
is rather dubious considering the important uncertainties 
of the model itself. 

4-2. Influence of the secondary irradiation 

Irradiation of the secondary enhances mass transfer. 



Hameury et al. (1997) showed that if one assumes that 
the effect of irradiation is given by equation (^ , outbursts 
having the general characteristics of superoutbursts (long 
durations, flat top or exponential decay with an abrupt 
cut-off) are expected. This model was, however, applied to 
a case in which the quiescent mass transfer was low enough 
for the disc to be stable on the cool branch; the instabil- 
ity was triggered by an external perturbation of the mass 
transfer from the secondary. This was required to explain 
the very long recurrence times of systems such as WZ Sge 
when standard values of a are assumed. Ma rgina lly unsta- 
ble mass transfer rates as in Warner et al. (1996) can give 
similar recurrence times but also in this case the amount 
of mass accreted during the superoutburst requires a sub- 
stantial enhancement of mass transfer. 

Figure ^ shows the light curve obtained when one in- 
cludes the secondary irradiation in the model. We neglect 
here disc irradiation and the disc is not truncated. We 
have taken 7 = 0.5, and all other parameters are as in 
Fig. (|ll), for a 1 Mq primary (i.e. Mq = 3 x 10^^ g s^^). 
The light curve is similar to those observed in SU UMa 
systems; it shows several normal outbursts separated by 
a large one which is sustained by enhanced mass transfer 
from the secondary. 

Large outbursts occur when the surface density at the 
outer edge of the disc is large enough that a cooling wave 
does not start immediately after the heating wave has ar- 
rived; equivalently, the disc mass must be larger than some 
critical value, and one therefore expects that the recur- 
rence time of such large outbursts Tg varies roughly as 
Mo . This, however, requires that the disc mass keeps 
increasing despite the presence of small outbursts, which 
means that Mq must be large enough to refill the disc 
with more mass than is lost during such outbursts. For 
Mo = 3 X 10"'^^ g s~^, which is more appropriate for short 
period systems, one does not get long outbursts, at least 
for the value of a considered here. Short outbursts are (as 
all our outbursts) of the inside-out type, and their recur- 
rence time Tn is the viscous time, and therefore do not 
depend on the mass transfer rate. 

The correlati on betw een (x T°-^ found in SU UMa 
systems (Warner 1995^) is interpreted, in the framework 
of the tidal-thermal instability (TTI) model, as result- 
ing from Ts varying as M~^, as in our case, and from 
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Fig. 3. Effect of secondary illumination on the predicted 
light curve; the parameters are those of Fig. |l|, with Mi ~ 
0.6 M0 and 7 = 0.5. The upper panel shows the outer disc 
radius, the middle panel shows the accretion rate onto the 
white dwarf (solid line) and the mass transfer rate from 
the secondary (dashed line); the lower panel shows the 
visual magnitude of the disc. 



Tn oc for outside-in outbursts (Osaki 1995a). One 

must however be careful with such a simple interpreta- 
tion. This explanation is valid only if M is the only pa- 
rameter determining both and Tg] this is clearly not 
the quantities such as the viscosity in quiescence 

(that may vary by orders of magnitude from WZ Sge type 
systems to "normal" SU UMa's), the disc radius (to the 
power 5.6), and the orbital period enter together with M 
in expressions for Tn and Tg (Osaki 1995a). It must also 
be stressed out that, for the low mass transfer rates of SU 
UMa stars, outside-in outbursts are not a natural outcome 
of the models and are produced by lowering the value of 
Qfc or by making it an appropriate function of radius. 

The time evolution of the outer disc radius is differ- 
ent from the predictions of the TTI model in several re- 
spects: (i) Tout varies during a normal outburst, whereas 
in the TTI model, rout remains roughly constant during 
an outburst after an initial increase during the rise (ii) 
''out oscillates at the beginning of a large outburst, (iii) 
the disk extends to a larger radius during a large outburst 



than during a short outburst, allowing for the possibility 
of the development of superhumps if the radius can reach 
the 3:1 resonance radius, whereas in the TTI model the 
disk size varies by 30% during a superout burst, with an 
average that is smaller than the average size during the 
previous normal outburst, and (iv) Tout remains approx- 
imately constant during a superoutburst, showing only a 
slow decline. 

Our results are similar to those of Smak ( 1991b| ), who 
considered the effect of enhanced mass transfer during 
superoutbursts, and concluded that observed disc radius 
variations in Z Cha and the length of the cycle in VW 
Hydri appeared to support the enhanced mass transfer 
model. The main difference with our work comes from 
the approximations describing the effect of illumination: 
whereas we assume a dependence between Mtr and Mace 
given by Eq. ^, in Smak's model Mtv is increased by about 
one order of magnitude after the maximum of a normal 
outburst, during a fixed period. As a consequence, the 
disc radius we obtain at the end of a large outburst is 
much smaller than in Smak's model. In our case, when 
the cooling wave starts propagating, the accretion rate 
onto the white dwarf, and hence the mass transfer from 
the secondary, is unaffected, and the disc contracts as in 
the unilluminated case, whereas in Smak's model, the disk 
expands rapidly when mass transfer is reduced by a fac- 
tor 10; the surface density then drops at the outer edge 
below the critical value, and a cooling wave starts in quite 
a large disc. 



Ichikawa et al. (1993) also considered a mass trans- 
fer outburst as the source of superoutbursts. They com- 
pared the resulting lightcurves with those produced by 
the tidal-thermal model and concluded in this last model 
gives a much better representation of observed properties 
of SU UMa's system. One should stress, however, that in 



Ichikawa et al. (1993) superoutbursts are triggered by a 



mass transfer 'instability'. In our case, superoutbursts are 
triggered by the usual thermal-viscous instability and it 
is only the subsequent evolution of the outburst which is 
modified by an enhanced mass transfer. In Ichikawa et al. 
(1993), the mass transfer Mtr is increased by a factor 100, 
whereas in our case the peak mass transfer rate from the 
secondary is 1.2 x 10^^ g s~^, i.e. increased only by a fac- 
tor 4. It is much smaller than the maximum possible rate 
from an irradiated low- mass star (see Section 2.1). Note 
that the average mass transfer rate, as given by Eq. ^ 
is larger than Mq, and thus larger than for Fig. |l|; it is in 
this case 4.8 x 10^^ g s^^ 

Observations are, for the moment, not of much help 
in deciding which of the models is right. There are good 
reasons to believe that a tidal instability is required to ac- 
count for the superhump phenomena. What is not known, 
however, is (i) the increase in the tidal torque resulting 
from this instability, and (ii) the radius at which the insta- 
bility stops. There is also good evidence for an enhanced 
mass transfer, due to irradiation, during outbursts but a 
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rehable description of this effect is missing. In both mod- 
els, however, the observed correlation with superout burst 
and normal outburst frequency can be obtained only by 
playing with the viscosity prescription which, of course, is 
not very satisfactory. 

The quiescent luminosities and accretion rates are al- 
most identical in the standard and enhanced mass transfer 
cases, as expected, since after the passage of the cooling 
wave, the disc has essentially forgotten its initial condi- 
tions; differences arise only from the different disk sizes 
which are smaller in the illuminated case because of the 
large mass transfer increase. 

4-. 3. Parameter dependence 
4.3.1. Mass transfer rate 

In this section we shall consider the combined effect of 
irradiation of both the disc and the secondary and shall 
determine how the resulting lightcurves depend on the 
value of rate AIq at which mass is transferred from an 
unilluminated secondary. 

Figure |4| shows various light curves obtained by varying 
Mq. We considered a system containing a 1 Mq primary, 
whose radius is 5 x 10^ cm and surface temperature 35,000 
K in quiescence. We have taken 7 = 0.5, and Mq ranges 
from 10"'^^ to 7 X 10^^ g s~^; this corresponds to average 
transfer rates in the range 1.5 - 8.6 xlO^^ g s^^. The inner 
disc radius is equal to the white dwarf radius. The white 
dwarf temperature has been chosen in the upper range of 
observed values in order to emphasize illumination effects. 
The effect is quite dramatic. Light-curves corresponding 
to cases (a) and (b) do not seem to be observed (as men- 
tioned in the previous section the mass transfer rate is too 
low for long outburts to be present). It might be that the 
corresponding systems exist, but have not yet been dis- 
covered because they are intrinsically rare - the parame- 
ters of Fig. ^ are at the upper range of allowed values, or 
that some of the curves we obtain are artifacts due to our 
oversimplified treatment of the secondary response to illu- 
mination. Lightcurves (c) and (d), however, compare very 
well with those of systems having very short supercycles 
such as RZ LMi. In our model they are obtained for high 
mass transfer rates, which is natural since these systems 
spend most of their time in the high state, whereas in the 
tidal-thermal instability such lightcurves require an ad 
hoc reduction of the parameter describing the tidal inter- 
action. The agreement with RZ LMi can be improved; in 
particular a reduction of the duration of the long outburst 
will be obtained by decreasing 7. 

One should note, however, that we have assumed that 
Mtr responds immediately to changes in Mace, whereas 
one could argue that there is a delay of the order of one 
or two days between illumination and the increase of mass 
transfer as discussed above. The introduction of such a de- 
lay makes the occurrence of long outbursts more difficult, 




50 100 
time (days) 



150 



Fig. 4. Visual magnitudes obtained in the case of a system 
containing a 1 M0 primary, whose radius is 5 x 10^ cm 
and surface temperature 35,000 K in quiescence. 7 = 0.5 
for aU four panels, and Mq = lO^^, 3 x lO^*', 4 x 10l^ and 
7 X 10^^ g s^^ from top to bottom. Note that the average 
mass transfer from the secondary is larger than Mq. 



as these require a near balance between mass transfer and 
accretion onto the white dwarf that must be established 
within a short outburst. We checked that if we use in Eq. 
(^ the average of the mass accretion rate over the past 
2.5 days (the duration of short outbursts), long outbursts 
are suppressed. 

Finally, it is worth pointing out that small outbursts in 
Fig. ^ are intrinsically different from those obtained when 
the disc illumination is not taken into account. Whereas 
in a non-irradiated disc small amplitude outbursts appear 
when the heating front cannot bring the whole disc into a 
hot state, here the disc never returns to quiescence in its 
inner parts, so the cooling wave is reflected into a heat- 
ing wave when it gets close to the stable hot inner part 
of the disc. The amplitude of these reflares grows as a 
consequence of the enhanced mass transfer during maxi- 
mum, until the disc mass has grown up to a point where 
a self-sustained long outburst is possible. 
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0.02; all other 



Fig. 5. Light curves obtained with acoid 
parameters are the same as for Fig. |b. 

4.3.2. Viscosity 

The reflares properties also depend on the ratio 
Q^hot/c^coid : the smaller this ratio, the more important the 



reflares (see Menou et al. 1999| for a discussion of this effect 
in the context of X-ray transients). This is simply due to 
the fact that, the lower this ratio, the larger S/Smax after 
the passage of a cooling front, E^ax being the maximum 
surface density on the cold stable branch; in the limiting 
case tthot = ttcoid there are no outbursts (Smak 1984| ), but 
a heating/cooling wave that propagates back and forth. 
Fig. H shows the effect of changing acoid to a smaller (by 
a factor 2) value. Successive reflares no longer reach the 
outer edge of the disc, and their amplitude therefore de- 
creases from one mini-outburst to the next one. This ac- 
counts for the presence of flat top outbursts which were 
absent in Fig. ^b. The lightcurves are similar for all mass 
transfer rates, showing the pattern of Fig ^ with longer re- 
currence times for smaller Aftr. The only exception is for 
7 X 10"'^^ g s~^, which is close to stability, and for which 
the main outbursts are of the outside-in type. 



4.3.3. Outer disc radius 

Since small discs favour large reflares, it is not surpris- 
ing that when one considers discs with average rout = 
1.3 X 10^° cm, and one takes acoid = 0.02 and ahot = 0.2, 
one obtains a combination of the lightcurves shown in the 
two previous sections; Fig ^ is a good example of this. It 
is worth noting that such a light curve is reminiscent of 
that of EG Cnc, even though the timescales are not quite 
the same. We do obtain the right pattern for the reflares, 
but we do not reproduce the very long superoutburst of 
EG Cnc (100 days), that would require 7 to be very close 
to unity, meaning that the linear approximation in Eq. ^ 
is invalid. For WZ Sge, one already had to assume a rel- 
atively large value of 7 (0.87) in order to reproduce the 
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Fig. 6. Example of a lightcurve obtained for a small disk; 
we consider a 1 M0 primary with surface temperature 
35,000 K, ttcoid = 0.02, ahot = 0.2, 7 = 0.5, and Mq = lO^^ 
g s"^ 



observed 25 days duration; since the outburst duration 
varies as l/log(7) (Hameury et al. 1997 ), we would need 
7 = 0.97 to obtain 100 days. Another difference with EG 
Cnc is the amplitude of the miniout bursts: the observed 
ones have approximately the same amplitude, whereas we 
get two identical miniout bursts, the others being of de- 
creasing amplitude. We have not been able to reproduce 
this behaviour with our parameterization; a possible solu- 
tion is to introduce a time dependent temperature of the 
white dwarf. This is expected, because the superoutburst 
lasted long enough to heat up the surface of the white 
dwarf that will then cool. 

If the rebrightenings of EG Cnc are indeed due to il- 
lumination effects, this implies that acoid cannot be small 
as we do not obtain refla res wh en Ocoid is significantly less 
than 0.01; Osaki et al. ( 1997 ) reached the same conclu- 
sion, but on different grounds; they assumed that acoid 
was increased to 0.1 during the superoutburst, remained 
high for 2 months, and then decreased back to small val- 
ues (0.001), and had therefore to set acoid to be an explicit 
function of time. 



4.3.4. Inner disc radius 

Finally we consider the effect of removing the inner disc 
regions, keeping both the disc and secondary irradiated. 

Apart from increasing the delay between the onset of 
an outburst in the disc and accretion onto the white dwarf, 
a large inner disc radius has a stabilizing effect on the disc 
itself, by preventing inside-out outbursts. This effect is 
quite noticeable in the case where the white dwarf surface 
temperature is high: if the inner disc radius rin is large 
enough, the unstable transition region between the stable 
region heated above hydrogen ionization temperature and 
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Fig. 7. Lightcurve obtained when the inner disc is trun- 
cated at a radius given by Eq. (|^), with a magnetic mo- 
ment of 2 X 10^° G cm^, the other parameters being that 
of Fig. |c. 

the cooler, quiescent external part does not exist. This 
suppresses the bounces after a longer outburst, as can be 
seen in figure 

For a given mass transfer rate, there is a critical value 
of the inner radius above which the disc is stable; when 
Tin approaches this value, the recurrence rate goes to in- 
finity. Arbitrarily long recurrence rates could therefore be 
expected, but only at the expense of very fine tuning; in 
the normal case where the mass accretion rate onto the 
white dwarf is negligible in quiescence as compared with 
the mass transfer rate from the secondary, the reasoning 
used by (Smak 1993) applies. The recurrence time irec is 
equal to: 

AM 



crit 



(8) 



Aftr Mtr 

where / is the ratio of the amount of mass transferred 
during an outburst AM and the maximum possible disc 
mass Merit, obtained assuming that the surface density is 
everywhere the critical surface density. Numerical models 
show that the surface density is not very far from its criti- 
cal value even at large radii, and that the amount of mass 
transferred during a normal outburst is typically 10% of 
the total disc mass. Therefore, / is not a very small pa- 
rameter that could freely vary, and large chang 
cannot result from variations in alone. A similar sit- 
uation is encountered in the case of soft X-ray transients 
(Menou et al. |1999[ ). 



5. Conclusions 

We have shown that many types of light curves can be 
produced by numerical models that include the illumina- 
tion of both the secondary and the accretion disc, thereby 
explaining a great variety of observed light curves. These 



effects account for phenomena such as post-outburst re- 
brightening (e.g. EG Cnc), long outbursts (U Gem for 
example), or SU UMa systems with extremely short su- 
percycles. In order to explore further these possibilities, 
one would need to determine from observations the mass 
transfer rate from the secondary as a function of the mass 
accretion rate onto the white dwarf, with a better accu- 
racy than it is available now. 

Despite the fact that our approximations are very 
crude, in particular the one concerning the response of 
the secondary to illumination, we can nevertheless draw a 
number of conclusions. First, the illumination of the disc is 
important only if the white dwarf is relatively massive, so 
that it can have a high temperature without contributing 
too much to the light emitted by the system in outburst, 
and that the efficiency of accretion is high. Rebrightenings 
also require a not to be too low in quiescence. 

The fact that we can reproduce an alternance of nor- 
mal and long outbursts when the illumination of the sec- 
ondary is included does not of course imply that the 
thermal-tidal instability model for SU UMa is incorrect; 
a tidal instability is most probably required to account 
for the superhump phenomenon. The question of the pre- 
cise role of this instability is, however, still open and our 
results raise some doubt on the validity of the parame- 
ters derived when fitting the observations, in particular 
for systems having a very short supercycle. 

This also means that the determination of the viscosity 
from the modeling of light curves is a far more difficult 
task than previously estimated. We obviously need some 
progress in the determination of the tidal torque; we also 
need to know how the secondary responds to illumination. 
We finally should include 2D effects in our models. First 
because the orbits in the outer disc are far from being 
circular, and second because the presence of a hot spot 
whose temperature can be of order of 10,000 K could in 
principle significantly alter the stability properties of the 
outer disc. 
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